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ABSTRACT 



A high fraction of the Low Mass X-ray Binaries (LMXBs) in early-type galax- 
ies are associated with globular clusters (GCs). Here, we discuss the correlations 
between LMXBs and GCs in a sample of four early-type galaxies with X-ray 
source lists determined from Chandra observations. There is some evidence that 
the fraction of LMXBs associated with globular clusters (/x-gc) increases along 
the Hubble sequence from spiral bulges (or spheroids) to SOs to Es to cDs. On the 
other hand, the fraction of globular clusters which contain X-ray sources appears 
to be roughly constant at /gc-x ~ 4%. There is a strong tendency for the X-ray 
sources to be associated with the optically more luminous GCs. However, this 
correlation is consistent with a constant probability of finding a LMXB per unit 
optical luminosity; that is, it seems to result primarily from the larger number 
of stars in optically luminous GCs. The probability of finding a bright LMXB 
per unit optical luminosity in the GCs is about L5 x 10~^ LMXBs per Lqj for 
Lx > 1 X 10^^ erg s"^ (0.3-10 keV) and rises to about 2.0 x 10"^ LMXBs per 
Lqj at lower X-ray luminosities, Lx ^ 3 x 10^'' erg s^^. This frequency appears 
to be roughly constant for different galaxies, including the bulges of the Milky 
Way and M31. There is a tendency for the X-ray sources to be found preferen- 
tially in redder GCs, which is independent of optical luminosity correlation. This 
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seems to indicate that the evolution of X-ray binaries in a GC is affected either 
by the metalUcity or age of the GC, with younger and/or more metal rich GCs 
having more LMXBs. There is no strong difference in the X-ray luminosities of 
GC and non-GC LMXBs. There is a weak tendency for the brightest LMXBs, 
whose luminosities exceed the Eddington luminosity for a 1.4 Mq neutron star, 
to avoid GCs. That may indicate that black hole X-ray binaries are somewhat 
less hkely to be found in GCs, as seems to be true in our Galaxy. On the other 
hand, there are some luminous LMXBs associated with GCs. There is no clear 
evidence that the X-ray spectra or variability of GC and non-GC X-ray sources 
differ. We also find no evidence for a difference in the spatial distribution of 
GC and non-GC LMXBs. Many of these results arc similar to those found in 
NGC 1399 and NGC 4472 by Angelini et al. and Kundu et al., respectively. 

Subject headings: binaries: close — galaxies: elliptical and lenticular — galaxies: 
star clusters — globular clusters: general — X-rays: binaries — X-rays: galaxies 



1. Introduction 

X-ray observations since the time of the Einstein Observatory have shown that early- 
type galaxies arc often luminous X-ray sources (e.g., Forman, Jones, & Tucker 1985). For 
the X-ray-luminous early-type galaxies (defined as those having a relatively high ratio of 
X-ray to optical luminosity Lx/Lb), it is clear that the bulk of the X-ray luminosity is from 
hot (~10^ K) interstellar gas. On the other hand, recent Chandra observations have resolved 
most of the X-ray emission in the X-ray-faint early-type galaxies into individual point-like 
sources (e.g., Sarazin, Irwin, & Bregman 2000, 2001). Given their properties and the stellar 
populations in these galaxies, these X-ray sources are assumed to be Low Mass X-ray Binaries 
(LMXBs). LMXBs are also observed in Chandra observations of X-ray- luminous early-type 
galaxies, and often dominate the X-ray emission at hard X-ray energies (e.g., Kraft et al. 
2000; Angelini, Loewenstein, & Mushotzky 2001). 

Chandra observations of early-type galaxies have shown that a significant fraction (~20- 
70%) of the LMXBs arc associated with globular clusters in the host galaxies (Sarazin et al. 
2000, 2001; Angelini ct al. 2001; Kundu et al. 2002). The fraction of LMXBs located in GCs 
is much higher than the fraction of optical light, which indicates that nondegenerate stars 
in GCs are much more likely (by a factor of ~300) to be donor stars in X-ray binaries than 
field stars. As has been known for a number of years, a similar result applies to our own 
Galaxy and to the bulge of M31 (Katz 1975; Hertz & Grindlay 1983; White et al. 1995). 
This is generally beheved to result from stellar dynamical interactions in globular clusters. 
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which can produce compact binary systems (Clark 1975; Fabian, Pringle, & Rees 1975; Hills 
1976). 

X-ray observations with ASCA indicate that the total luminosity of LMXBs in early- 
type galaxies correlates better with the number of GCs than with the optical luminosity of 
the galaxy (White, Sarazin, & Kulkarni 2002). The correlation is consistent with the total 
X-ray luminosity of LMXBs being proportional to the number of GCs. This is somewhat 
surprising, as a nontrivial fraction (~ 50%) of the LMXBs in most of the early-type galaxies 
observed so far with Chandra are not identified with GCs. This suggests that most (perhaps 
all?) of the LMXBs in carly-typc galaxies were made in GCs (Grindlay 1984; Sarazin et al. 
2001; White et al. 2002). The field LMXBs might have been ejected from globular clusters by 
kick velocities resulting from supernovae, by stellar dynamical processes, or by the dissolution 
of the globular due to tidal effects. If this is true, then the field LMXBs are a useful record 
of the total globular cluster populations and of the history of globular cluster destruction in 
galaxies. 

Alternatively, if field LMXBs are not generally made in GCs, then one might expect 
field LMXBs and GC LMXBs to have different properties. Presumably, the stars in field 
LMXBs are formed in binaries, and evolve independently (i.e., they are not affected by their 
environment). As evidenced by the high X-ray to optical ratio of GCs, LMXBs in GCs 
are strongly affected by dynamical interactions within the GC. For example, it has been 
argued that it is difficult for GCs to retain LMXBs with massive black holes (BHs), which 
may be ejected by dynamical interactions (e.g., Portegies Zwart & McMillan 2000), which 
is consistent with the fact that there are presently no BH LMXBs identified with GCs in 
our own Galaxy. Thus, one might expect very luminous LMXBs to avoid GCs in early-type 
galaxies. 

One might also expect LMXBs in different GCs to be affected by differences in their GC 
environment. This might occur because of differences in the ages or metallicities of different 
GCs, which might affect the evolution of individual LMXBs. Alternatively, differences in 
stellar densities would change the rates of stellar dynamical interactions. 

In this paper, we give a preliminary analysis of the correlations of GC LMXBs and 
non-GC LMXBs, and of the differences between X-ray GCs and non- X-ray GCs in early- 
type galaxies. This analysis should be viewed as preliminary because of the small size and 
incompleteness of the existing samples. At present, only a relatively small number of early- 
type galaxies have been observed with Chandra. The spatial resolution of Chandra is needed 
to separate LMXBs at the distance of the Virgo cluster, where most of the nearby early-type 
galaxies are located. Moreover, most of the galaxies observed with Chandra so far are X-ray 
bright galaxies. In these galaxies, bright diffuse emission from hot gas makes it difficult to 
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detect LMXBs near the center of the galaxy (e.g., Angehni et al. 2001). The optical samples 
of GCs also are incomplete. Ground-based observations have difficulty in detecting GCs near 
the centers of early-type galaxies, where they are blended into the diffuse optical emission 
from the galaxy. Thus, ground-based GC samples tend to be restricted to the outer parts 
of the galaxies, which is not where most of the LMXBs are located. Also, ground-based 
observations do not generally resolve GCs at the distance of the Virgo cluster. Thus, there is 
a significant chance that ground-based GCs detected in the outer parts of early-type galaxies 
may actually be foreground or background objects. GC observations with the Hubble Space 
Telescope (HST) are greatly preferable, as the resolution of HST allows GCs to be detected 
in the inner regions of galaxies, and slightly resolves them. However, some of the galaxies 
with the best X-ray data have no HST GC observations (e.g., NGC 4697). Also, the field of 
view of a single Wide Field Planetary Camera 2 (WFPC2) observation is smaller than the 
angular size of early-type galaxies at the distance of the Virgo cluster. Multiple pointings 
with the WFPC2 (e.g., Kundu ct al. 2002) or with the recently installed Advanced Camera 
for Surveys (ACS) may improve this situation. 

Because the GC observations of most galaxies do not cover the entire galaxy, we will 
restrict our analysis to the LMXBs located in the regions covered by the GC sample. We will 
also restrict our discussion to the relative properties of GC LMXBs and non-GC LMXBs, 
and of X-ray GCs and non- X-ray GCs. Hopefully, comparing the relative properties of these 
samples within the spatial regions they share should reduce most of the selection effects in 
these incomplete samples, particularly as the LMXBs were detected solely by their X-ray 
properties and the GC were detected solely by their optical properties. All other things 
being equal, there is no reason why a LMXB in a GC is any more or less detectable in 
the X-ray than a field LMXB. Because of the irregular spatial regions covered by the GC 
samples, we will concentrate on the intrinsic properties of the sources rather than their 
spatial distribution. 

Our galaxy sample. X-ray source lists, and the detection limits for the LMXBs are 
discussed in § 2. The statistics of occurrence of LMXBs in the galaxies are given in § 3. 
The relative properties of GC and non-GC LMXBs are analyzed in § 4, while the relative 
characteristics of X-ray and non-X-ray GCs are determined in § 5. In § 6, we briefiy discuss 
the relative spatial distributions of the samples. We discuss our results and summarize our 
conclusions in § 7. 
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2. Sample 

Table 1 lists some of the properties of the early-type galaxies in our sample. The first 
five columns give the name, Hubble type, Galactic absorbing column Nh (Dickey & Lockman 
1990), effective radius rgg, and distance D. For consistency, all of the distances for these 
galaxies are taken from Tonry et al. (2001), based on the method of surface brightness 
fluctuations. For these galaxies, the statistical error in their distance modulus is generally 
about 0.17, leading to an error in the distance of 8%. There also is a systematic error (which 
would affect all of the galaxies by the same factor) of about 8% in the distance. 

The last 5 columns in Table 1 give the useful exposure time in the Chandra X-ray image, 
the minimum detectable source count rate, and the minimum detectable X-ray luminosity in 
the 0.3-10 keV band, and the references for the original papers giving the detections of the 
LMXBs and the GCs. The criteria for the detection of the X-ray sources are discussed in 
detail in our original papers on the X-ray observations. These papers also list the properties, 
including the count rates, for each of the sources. In every case, a wavelet detection algorithm 
(CIAO WAVDETECT^ program), was used to detect the sources, with a significance threshold 
of 10^®. This implies that <1 false source (due to a statistical fiuctuation in the background) 
would be detected in the entire S3 image. We further restricted the sources by requiring 
sufficient source counts to determine the source flux at the > 3(7 level. In general, this 
implies that the minimum detectable sources had >10 net counts. 

For the Virgo cluster galaxies, some of the X-ray sources have been identified with outer 
globular clusters from the lists given by Hanes (1977). In general, less information is available 
for these outer globulars. Also, there is a higher chance that the globulars and X-ray sources 
at large distances from the centers of the galaxies are actually foreground or background 
objects. Thus, we haven't included these globulars and the associated X-ray sources in our 
sample. 

Table 2 hsts some of properties of the X-ray sources in the region covered by the GC 
surveys. More detailed information is given in the original references on the X-ray detections, 
which are listed in Table 1. Table 2 gives the galaxy, the source number from the tables in 
the original references, the net X-ray count rate and error in the 0.3-10 keV band, the X-ray 
luminosity Lx in the same band (§ 5.1 below), the X-ray hardness ratios H21^ and iJ31*^ 
(§ 5.2 below), an indication of whether the source was found to be variable (§ 5.3 below), a 
check if the source was identified with a GC, the I-band absolute magnitude Mj of the GC 
(§ 4.1 below), and the V — I color of the GC (§ 4.2 below). The last column has a check if 



^See http://asc.harvard.edu/ciao/. 



-6- 



the source belongs to an X-ray luminosity limited complete sample which is defined in § 2.5 
below. 

Figure 1 shows the regions surveyed to detect GCs in these galaxies. These are compared 
to the circle with a radius of one effective radius. 

2.1. NGC 1553 

The globular cluster identifications in NGC 1553 come from a survey of globulars in 
SO galaxies with the HST Wide Field / Planetary Camera 2 (WFPC2; Kundu & Whitmore 
2001b). This GC sample included the nucleus of the galaxy and the inner parts of the galaxy 
to the southeast of the nucleus, out to a radius of about 2' (Fig. 1). 

The X-ray sources were published in Blanton et al. (2001). We adjusted the values of 
the X-ray luminosities to be consistent with the new distance to this galaxy (Tonry et al. 
2001). We excluded the central X-ray source (Src. 1, CXOU J041610.5-554646) because 
it is the central active galactic nucleus (AGN). Srcs. 2 (CXOU J041609.9-554646) and 6 
(CXOU J041611. 7— 554652) were excluded because they arc located near the center of the 
galaxy and 1-2" from several globular cluster positions. This is too large a separation for 
a definite identification in this crowded region, so we were uncertain whether to associate 
these sources with the GCs or not. 

The globular cluster identifications are discussed in Blanton et al. (2001). We used 
the position of the central AGN, which is present in both the HST optical and Chandra 
X-ray images, to register the two images, assuming a simple shift. This required an offset of 
1'.'9. [We also registered the images individually by comparison to accurate optical positions 
from the U.S. Naval Observatory A2.0 optical catalog (USN0A2 Monet D., et al. 1998) and 
the Two Micron AU Sky Survey (2MASS Cutri et al. 2001). This indicated that both the 
Chandra and iJST positions were off by ~1'.'5, and that the relative offset was the combination 
of these systematic errors.] Three potential GCs were found to agree in position with LMXBs 
to within 0'.'5. Based on simulations of the distributions of GCs and LMXBs, the expected 
number of random associations would be about 0.05. However, one of the potential GCs had 
a very blue color, and was dropped from the sample using the criteria for GC identifications 
given in Kundu & Whitmore (2001b). 



2.2. NGC 4365 

Sivakoff, Sarazin, & Irwin (2003) have determined the X-ray source population in 
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NGC 4365. The globular cluster population is from an HST WFPC2 survey of elliptical 
galaxies by Kundu & Whitmore (2001a). This GC sample included the nucleus of the galaxy 
and the inner parts of the galaxy to the northwest of the nucleus, out to a radius of about 
2' (Fig. 1). We found an offset of about 1'.'3 between the HST and Chandra positions. After 
removing this offset, 18 X-ray sources were found to lie within 1'.'02 of a GC. Simulations of 
the distributions of GCs and LMXBs indicate that about 0.6 associations would be expected 
at random. One additional association was found between a LMXB and a GC in the short 
list of clusters with spectroscopy recently published by Larsen et al. (2003). However, it 
is outside of the region covered by the GC survey in Kundu & Whitmore (2001a), and we 
haven't included this source in the present study. 



2.3. NGC 4649 

The X-ray sources in NGC 4649 (M60) have been published in Randall, Sarazin, & Irwin 
(2003). As with NGC 4365, the globular cluster sample here is from Kundu & Whitmore 
(2001a). This sample included the nucleus of the galaxy and the inner parts of the galaxy 
to the northwest of the nucleus, out to a radius of about 2' (Fig. 1). At the northwest 
corner of the WFPC2 field, there is a high density of globular clusters which are probably 
associated with the companion Sc galaxy NGC 4647. Because we wish to study the LMXB 
population in early-type galaxies, we have excluded this region from our sample. We found 
an offset of O'.'Q between the HST and Chandra positions. After correcting for this offset, 22 
GCs associated with NGC 4649 were found to coincide with X-ray sources to better than 
0'.'87. Based on simulations of the distributions of GCs and LMXBs, the expected number of 
random associations would be about 2.96, so a few of these identification may be spurious. 
NGC 4649 is an X-ray bright eUiptical galaxy (the only one in our sample). Unfortunately, 
the bright X-ray emission by interstellar gas makes it difficult to detect faint sources near 
the center of the galaxy, and the LMXB source list is incomplete within approximately 70" 
of the center of NGC 4649. This X-ray incompleteness may reduce the detected fraction 
of GCs containing X-ray sources, /gc-x- Ideally, one would exclude this region from our 
sample, but this would leave almost no area of the galaxy which was on the WFPC2 image, 
beyond 70", and not confused with the companion NGC 4647. Since our primary aim is a 
comparative study of globular cluster and non-globular cluster LMXBs, we have included 
all of the covered area of the galaxy. We exclude X-ray Src. 1 (CXOU J124340.0+113311), 
which is an extended feature due to structure in the diffuse gaseous emission. 
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2.4. NGC 4697 

The X-ray source population of NGC 4697 has been analyzed and listed in Sarazin et 
al. (2000, 2001). We adjusted the values of the X-ray luminosities to be consistent with the 
new distance to this galaxy (Tonry et al. 2001). The X-ray observation was done in Chandra 
Cycle- 1 at a time when significant absolute errors in positions occurred. Thus, the X-ray 
image was registered by comparing the positions of X-ray sources with optical identifications 
(not the GC discussed below) to optical positions from the USN0A2 catalog (Monet D., 
et al. 1998). This required an offset of about 2'.'6. Following this correction, the X-ray 
and optical positions agreed to better than 0'.'3. Subsequently, 2MASS (Cutri et al. 2001) 
positions became available for several X-ray sources with optical IDs, and these confirmed 
the offset and corrected accuracy of the X-ray positions. 

The globular cluster sample for NCC 4697 comes from an unpublished, ground-based 
ccd survey very kindly provided by J. Kavelaars (2000, private communication). These 
globulars are all located in a circular annular region extending from about 1'5 to 2'.5 from 
the center of NGC 4697 (Fig. 1). Seven of the X-ray sources were found to lie within 
0'.'78 of X-ray sources. From simulations based on the distribution of LMXBs and GCs, we 
would expect 0.18 associations to occur at random. Although the GCs are not resolved in 
these ground-based observations, this sample doesn't extend to a large enough projected 
distance that an important fraction of the GCs are hkely to actually be background AGN 
or foreground stars. It is difficult to detect globular clusters in the central regions of early- 
type galaxies with ground-based observations because of the effects of seeing and the high 
surface brightness of the galaxy. The available data on NGC 4697 do not include colors for 
the globular clusters, so these GCs will not be used for looking for trends with color. Also, 
the magnitudes for the NGC 4697 GCs were determined in the R-band, whereas the three 
other galaxies used the V and I bands. For comparison, we converted the R magnitudes to 
I magnitudes assuming that I — R — 0.49, which is appropriate for the average color of the 
GCs in the other galaxies. We also increased the errors to include a systematic error on this 
conversion. 



2.5. Combined Samples 

In this paper, we will discuss some of the properties of the LMXBs and associated GCs 
in the individual galaxies. However, our primary aim is to compare LMXBs and GCs in the 
sample formed by combining the galaxies. We will consider two combined samples. First, we 
will consider an X-ray luminosity hmited sample {"Lx complete") of sources whose X-ray 
luminosities exceed 1.23 x 10^^ erg s~^ in the 0.3-10 keV band, which is the highest detection 
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limit for the galaxies in the sample (Table 1). In order to improve the statistics, we will 
also consider the sample formed by combining all of the sources, although this sample is not 
complete in the X-ray luminosity. However, we will restrict our discussion to the relative 
properties of GC LMXBs and non-GC LMXBs, and of X-ray GCs and non-X-ray GCs. 
Hopefully, comparing the relative properties of these samples within the spatial regions they 
share should reduce most of the selection effects in this incomplete sample, particularly as 
the LMXBs were detected solely by their X-ray properties and the GC were detected solely 
by their optical properties. All other things being equal, there is no reason why a LMXB in 
a GC is any more or less detectable in the X-ray than a field LMXB. 



3. Statistics of Occurrence of LMXBs 

The numbers of LMXBs and GCs in our sample are listed in Table 3. The first three 
columns give the galaxy name, the Hubble type, and the number of GCs in the surveyed 
region A^gc- The next five columns give the number of LMXBs in the same region A^lmxb, 
the number of LMXBs coincident with GCs Aboth, the fraction of LMXBs identified with GCs 
fx-GC = (Aboth/ALMxe), the fraction of GCs associated with LMXBs /gc-x = (Aboth/A^Gc), 
and the number of LMXBs in GCs divided by their total optical luminosity. The LMXBs 
included in these columns are only those in the Lx complete sample {Lx > 1-23 x 10"^^ 
erg s~^ in the 0.3-10 keV band). The errors listed are l-o"; errors on fractions are derived 
assuming a binomial distribution. The next five columns repeat the statistics for the full 
sample of LMXBs. 

Note that the elliptical galaxies in our sample all have a similar fraction of LMXBs in 
GCs (/x-GC ~ 50%). The value for NGC 4697 in the complete sample is larger, but the 
errors are very big because of the small number of sources. In any case, the hypothesis that 
the number of LMXBs in GCs in NGC 4697 was drawn from the same distribution as that 
of the other ellipticals can only be rejected at the 57% confidence level. In the full sample, 
the fraction in NGC 4697 is in good agreement with those of the other eUipticals. A similar 
fraction was found for the X-ray bright eUiptical NGC 4472 by Kundu et al. (2002). 

Although the statistics are very poor, the fraction in the one SO galaxy in our sample 
(NGC 1553) is lower, fx-GC = 18%. We tested the hypothesis that the number of LMXBs 
in GCs in NGC 1553 was drawn from the same distribution as that for the elliptical galaxies, 
and found that it could be rejected at the 99% confidence level. The fraction in spiral bulges 
or spheroids appears to be even lower, /x-gc ~ 10% (e.g., Verbunt & van den Hcuvcl 1995). 
For example, in our Galaxy there are 14 LMXBs which have been associated with GCs 
(White et al. 1995; Harris 1996; Liu et al. 2001; White & Angehni 2001) out of a total of 
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about 150 LMXBs. This would suggest that /x-gc ~ 9%. However, most of the LMXBs 
seen in our Galaxy are much too faint to have been detected in distant galaxies. In the Milky 
Way, there are very few LMXBs with persistent X-ray luminosities as high as those for our 
Lx complete sample. So, we instead compare to the fraction seen in NGC 4697, where the 
luminosity limit is Lx > 2.7 x 10'^^ erg s~^. If we limit the sample to sources which are 
more luminous than this for periods of many hours (i.e., not the burst luminosity), then the 
fraction is very roughly /x-gc ~ 7%. At similar luminosity levels, the fraction of bulge or 
spheroid X-ray sources in M31 which are identified with GCs is about 20% (Primini, Forman, 
& Jones 1993; Supper et al. 1997; Barmby & Huchra 2001; Di Stefano et al. 2002). On the 
other hand, Angehni et al. (2001) report a fraction of ~70% in NGG 1399, the cD galaxy 
at the center of the Fornax cluster. Although the number of galaxies tested so far is small 
and the statistics in most of the galaxies are poor, all of this is consistent with an increase 
in /x-GC along the Hubble sequence from spiral bulges to SOs to Es to cDs. The specific 
frequency of globular clusters (the number per unit galaxy luminosity) also increases along 
the same sequence (e.g., Harris 1991). 

There is no clear trend in the fraction of GCs containing X-ray sources (/gc-x) with 
Hubble-type. All of our galaxies are roughly consistent with about 3.3% of the GCs hosting 
luminous LMXBs at the luminosity level of the Lx complete sample {Lx > 1.23 x 10'^^ erg 
s~^ in the 0.3-10 keV band). At lower luminosity levels, the fraction appears to increase to 
about 4%. Similar numbers are found for the spiral bulge or spheroid of our Galaxy and 
M31 (e.g., Harris 1996; Barmby & Huchra 2001; Liu et al. 2001; Di Stefano et al. 2002). If 
we compare to the fraction seen in NGG 4697, where the luminosity hmit is Lx > 2.7 x 10^^ 
erg s~^, we find that roughly /gc-x ~ 3%. Since many of the LMXBs are highly variable, 
and catalogs for our Galaxy include sources detected at different times, it is not clear how 
many of these sources could be detected in a single ^10 hour observation with Chandra if 
our Galaxy were at the distance of the Virgo cluster. In the bulge of M31, the fraction of 
GCs containing LMXBs with Lx > 2.7 x 10^^ erg s"^ is about 2% (e.g., Barmby & Huchra 
2001; Di Stefano et al. 2002). Similar fractions have been found for the eUiptical galaxy 
NGG 4472 (Kundu et al. 2002), and for the galaxy NGG 1399 at the center of the Fornax 
cluster (Angehni et al. 2001). In fact, Kundu et al. (2002) argued that /gc-x was constant 
at about 4% in a wide range of galaxies. This fraction may provide a constraint on the duty 
cycle of quiescence vs. X-ray brightness in the LMXBs in GCs. 

One concern with the fraction /gc-x is that it may be affected by incompleteness in the 
detection of faint GGs. Luminosity functions of globular clusters in early-type galaxies are 
often well fit by log-normal functions (e.g., Harris 1991); for such a luminosity distribution, 
faint GC may contribute to the number of GCs, but do not add significantly to the total 
luminosity. Below, we will find that LMXBs are preferentially associated with more luminous 
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GCs, but that the probabihty of finding a LMXB in a GC is roughly proportional to its 
optical luminosity. Thus, a more useful measure of the abundance of LMXBs in GCs may 
be the number per unit optical luminosity, which should be less sensitive to incompleteness 
or to the details of the luminosity function of GCs in each galaxy. This is given as the 
last numerical column in Table 1. For the entire sample of four galaxies, the total value is 
1.5 X 10"^ LMXBs per Lqj in the Lx complete sample and 2.0 x 10~^ LMXBs per Lqj in 
the full sample. To within the errors, all of the galaxies are consistent with this value. At 
luminosities Lx ^ 2.7 x 10^''' erg s~^ (2-10 keV), the corresponding value for our Galaxy is 
about 1.4 X 10-^ LMXBs per L©,/. 



4. GC Properties of LMXBs 
4.1. Luminosities 

Figure 2 shows histograms of the absolute I magnitude, Mj, of the total GC sample 
(upper histogram) and of the GCs containing LMXBs (shaded histogram) for the Lx com- 
plete sample (left) and the full source sample (right). The results for the individual galaxies 
NGC 4365 and NGC 4649 are shown in Figure 3. In NGC 1553 and NGC 4697, there are 
too few LMXBs associated with GCs to allow a useful comparison. All of the magnitudes 
were corrected for Galactic extinction. 

The LMXBs seem to be associated preferentially with the more optically luminous GCs. 
For example, the median values of Mj for non- X-ray GCs are —8.7 {Lx complete sample), 
-8.7 (full sample), -9.1 (NGC 4365). and -8.8 (NGC 4649), while the corresponding values 
for the X-ray GCs are -10.2 (Lx complete sample), -10.2 (fuU sample), -10.1 (NGC 4365), 
and —10.3 (NGC 4649). Using the Wilcoxon or equivalent Mann- Whitney rank-sum tests 
(Mann & Whitney 1947), the distribution of X-ray and non- X-ray luminosities are found 
to disagree at more than the 5a {Lx complete sample), 6a (full sample), Aa (NGC 4365), 
and 5a (NGC 4649), with the probability that they are drawn from the same distribution 
being kW^ {Lx complete sample), KlQ-^^ (full sample), <10-^ (NGC 4365), and <10-^ 
(NGC 4649). We also compared the distributions of X-ray and non-X-ray GC absolute 
magnitudes using the Kolmogorov-Smirnov (KS) two-sample test. The probability that 
they were drawn from the same distribution is less than 10~^ {Lx complete sample), 10~^ 
(full sample), lO'^ (NGC 4365), 10"^ (NGC 4649). Thus, the X-ray sources are strongly 
correlated with the brighter globular clusters. A similar result was found previously for the 
galaxy NGC 1399 by Angehni et al. (2001) and for NGC 4472 by Kundu et al. (2002). 

Of course, a correlation between optical luminosity and the probability of having an X- 
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ray source is not unexpected. LMXBs contain normal stars, and globular clusters which have 
higher luminosities have more stars as potential donors in LMXBs. Thus, it is interesting 
to test the hypothesis that the probability that a GC contains a LMXB is proportional to 
its optical luminosity. Figure 4 compares the cumulative probability distribution of LMXBs 
versus the cumulative distribution of the optical luminosity in GCs (for the Lx complete 
sample). Both were accumulated starting at small M/ (high optical luminosity) and going 
to higher Mj (low luminosity) . The optical luminosity is determined solely from the I-band 
absolute magnitude. The two cumulative distribution functions track one another fairly well. 
For example, half of the optical luminosity comes from GCs brighter than Mj — —10.1, 
while the medium absolute magnitude of GCs with LMXBs is —10.2. There is a slight 
tendency for the fainter GCs to have a higher probability to have an X-ray source per unit 
luminosity, but it is not statistically significant. However, the GC samples are incomplete 
at the faint end, and the completeness limit differs significantly from galaxy to galaxy in 
our sample. Although there is no obvious reason why this incompleteness would affect 
GCs with X-ray sources differently than those without LMXBs, this suggests caution in 
interpreting differences among faint GCs. We used a modified version of the KS two-sample 
test to compare the two distributions. The largest vertical difference in the cumulative 
distributions \s D — 0.115, and this difference is not significant given the number of LMXBs. 
Very similar results were found for the full sample (albeit more strongly) and for NGC 4365 
and NGC 4649 (but less strongly). 

Thus, the current data indicate that optically bright GCs are much more likely to contain 
LMXBs than faint GCs, but the distribution is consistent with a constant probability per 

unit optical luminosity. A similar result was found previously for the galaxy NGC 4472 by 
Kundu et al. (2002). As noted in § 3, the probability of having a LMXB per unit optical 
luminosity of the GC also appears to be constant from galaxy to galaxy, with a value of 
about 1.5 X 10"^ LMXBs per Lqj for > 1 x 10^^ erg s'^ (0.3-10 kcV). The value rises 
to about 2.0 X 10~^ LMXBs per Lqj at lower X-ray luminosities {Lx ^ 3 x 10^^ erg s"^). 

4.2. Optical Colors 

Figure 5 shows histograms of the V — I colors for the total GC sample (upper histogram) 
and for the GCs containing LMXBs (shaded histogram) for the Lx complete sample (left) 
and the full source sample (right). The colors are all corrected for Galactic reddening. 
Because these samples contain a mixture of three galaxies, the overall color distribution 
may be less obviously bimodal than that seen in some individual galaxies (e.g., Kundu et 
al. 2002). The results for the individual galaxies NGC 4365 and NGC 4649 are shown in 
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Figure 6. As has been noted previously (e.g., Kundu et al. 2002), NGC 4649 has a bimodal 

color distribution, whereas NGC 4365 was a broad color distribution which is not clearly 
bimodal. In NGC 1553, there are too few LMXBs associated with GCs to allow a useful 
comparison, and no colors were available to us for the GCs in NGC 4697. 

The LMXBs appear to be associated preferentially with the redder GCs (larger values 
of y — /), although the trend is not as strong as that for the absolute magnitude or optical 
luminosity (§ 4.1). The median colors of the non- X-ray GCs are V — I = 1.07 {Lx complete 
sample), V-I = 1.07 (full sample), V-I = 1.05 (NGC 4365), and V-I = 1.09 (NGC 4649), 
while the corresponding medians for the X-ray GCs are 1.12 {Lx complete sample), 1.14 (full 
sample), V -I = 1.09 (NGC 4365), and V - I = 1.21 (NGC 4649). Using the Wilcoxon or 
Mann- Whitney rank-sum tests, the probabilities that the distribution of X-ray and non-X-ray 
optical colors were drawn from the same distribution are 0.018 {Lx complete sample), 0.0019 
(full sample), 0.046 (NGC 4365), and 0.0069 (NGC 4649), which correspond to differences 
of 2.1a {Lx complete sample), 2.9a (fuU sample), 1.7a (NGC 4365), and 2.5a (NGC 4649). 
Using the KS test, the probabilities that the two color distributions where drawn from the 
same distribution are 0.030 {Lx complete sample), 0.0013 (full sample), 0.14 (NGC 4365), 
and 0.031 (NGC 4649). Thus, a significant trend for LMXBs to be associated with red 
GCs appears in the Lx complete sample, in the full sample, and in NGC 4649. A similar 
correlation with the optical color of GCs was found previously for the galaxy NGC 1399 by 
Angelini et al. (2001) and for NGC 4472 by Kundu et al. (2002). The trend is less obvious 
in NGC 4365. Recent optical/IR observations of GCs in NGC 4365 indicate that this galaxy 
has both old and intermediate age (~5 Gyr) red GCs (Puzia et al. 2002; Larsen et al. 2003), 
which may account for its broad but unimodal color distribution. The LMXBs seem to be 
preferentially associated with the younger red GCs (Sarazin et al. 2003). 

A correlation between LMXBs and GC color might result from the correlation between 
LMXBs and the optical luminosity of GCs, if the brighter globular clusters were preferen- 
tially redder. Figure 7 shows the scatter diagram between the optical colors and absolute 
magnitudes for the GCs, with those associated with LMXBs in the Lx complete sample indi- 
cated. Various correlation measures show no evidence for a significant statistical correlation 
between absolute magnitude and color. This test doesn't take into account the various selec- 
tion effects which went into constructing this sample; it only shows that there is no strong 
correlation between absolute magnitude and color within this sample. However, previous 
detailed optical studies, including selection effects, on the GCs in these and other early-type 
galaxies have found no significant correlation of color with optical luminosity (e.g., Kundu 
& Whitmore 2001a). 

The correlations between X-ray sources and GC color would seem to suggest that the 
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evolution of X-ray binaries in GCs is affected eitfier by the metallicity or age of the GC, with 
younger and/or more metal-rich GCs having more LMXBs. 

5. X-ray Properties of GCs 

5.1. X-ray Luminosities 

Figure 8 shows histograms of the numbers of X-ray sources as a function of their X-ray 
luminosities for all of the sources (upper histogram) and for the sources associated with GCs. 
This figure is for the full sample, but the lower limit luminosity for the Lx complete sample 
is also marked. Figure 9 shows the distribution of X-ray luminosities for the individual 
galaxies NGC 4365 and NGC 4649; NGC 1553 and NGC 4697 had too few GC-LMXBs to 
provide useful comparisons. The X-ray luminosities are in the band 0.3-10 keV and have 
been corrected for absorption. Although the statistical accuracy is hmited by the small 
number X-ray sources, at most luminosity levels there is no obvious difference in the two 
distributions. However, there may be a tendency for the most X-ray-luminous sources to 
avoid GCs. In the Lx complete sample, the rank sum test and KS test both find that the 
GC sources are fainter. The probability that the two populations are drawn from the same 
distribution can be rejected at the 97% level, corresponding to a difference of about 1.9(7. 

It is generally found that the luminosity functions of X-ray sources in early-type galaxies 
have a broken power-law form, with a break luminosity of ~ 3 x 10^^ erg s~^ (e.g., Sarazin 
et al. 2001). This break luminosity is close to the Eddington luminosity of a 1.4 Mq neutron 
star. As a result, Sarazin ct al. (2000, 2001) have argued that the sources with Lx > Lb 
predominantly contain black holes, while those below this luminosity may contain neutron 
stars. In the Lx complete sample. 19±6% of LMXBs in GC have luminosities greater than 
3 X 10^^ erg s~^, while the fraction for the non-GC sources is 44±8% {la error bars). Thus, 
the difference is somewhat significant. 

Our results are similar to those found for the X-ray bright elhptical galaxy NGC 4472 
by Kundu et al. (2002), who found no significant difference in the X-ray luminosities of GC 
and non-GC sources. On the other hand, Angehni et al. (2001) found that the brightest 
X-ray sources in NGC 1399 were mainly associated with GCs, which is the opposite of the 
weak tendency we find. 
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5.2. X-ray Heirdness Ratios 

We used hardness ratios or X-ray "colors" to crudely characterize the spectra of the 
sources. We defined two hardness ratios as 

«31 . (2) 

where S, M, and H are the total counts in the soft (0.3-1 keV), medium (1-2 keV), and 
hard (2-10 keV) bands, respectively. Although all of the galaxies lie in directions with low 
Galactic absorbing columns Nh (Table 1), there is a small variation in the absorption towards 
the different galaxies. Thus, we have corrected the hardness ratio to remove the effects of 
Galactic absorption. Let 5*0, Mq, and Hq be the soft, medium, and hard band counts, 
corrected for absorption. We define /21 = {S/ Sq)/{M/Mo), and /31 = {S/ Sq) / {H / Hq). 
Then, the absorption-corrected hardness ratios are 

^ (l + g21)/21-(l-g21) 

(1 + //21)/21 + (1 - i721) ' ^ ' 

rrn.O ^ {1 + H31) f31 - {1 - H31) 

{1 + H31)f31 + {1-H31) ■ ^ ^ 

We calculated the corrections assuming a single spectrum for all of the sources, which we 
took to be 7 keV thermal bremsstrahlung. The corrections are smaller than the errors in 
the hardness ratios for almost all of the sources, and the differences in the corrections from 
galaxy to galaxy are even smaller. 

In Figure 10, we show the absorption-corrected hardness ratios for the sources with more 
than 20 net counts in the Lx complete sample. The filled squares are the sources identified 
with GCs, while the open squares are non-GC sources. For comparison, the solid line shows 
the hardness ratios for power-law spectral models; the triangles indicate values of the power- 
law photon number index of F = (upper right) to 3.2 (lower left) in increments of 0.4. 
Although there are significant variations in spectra from source to source, on average the 
source spectrum can be fit by either thermal bremsstrahlung with a temperature of ~7 keV 
or a power-law with a photon index F ~ 1.5 plus Galactic absorption. The hardnesses ratios 
for the sources in the galaxies NGC 4365 and NGC 4649 are given separately in Figure 11. 
NGC 1553 and NGC 4697 had too few GC-LMXBs to provide useful comparisons. 

Histograms of the if21° and H31^ hardness ratios are given in Figure 12 for the Lx 
complete sample. Considering both Figures 10 and 12, no very significant differences in the 
hardness ratios are evident. The rank-sum and KS tests suggest that the H21^ colors of 
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the GC sources are slightly harder (larger if21°), but the difference is only significant at 
the 93% confidence level. The same tests indicate that the H31^ distributions are identical. 
Thus, there is no strong evidence that the spectra of the GC and non-GC sources differ in 
any systematic way. This agrees with the conclusions of Angelini et al. (2001) and Kundu 
et al. (2002) for NGC 1399 and NGC 4472, respectively 

There arc two X-ray sources (one GC and one non-GC source) with hardness ratios of 
(if21°, iJ31°) ~ (1,1)- It is likely that these X-ray sources are strongly self-absorbed. In 
some cases, such sources are found to be identified with background AGNs (e.g., Sarazin et 
al. 2000). However, both of these X-ray sources are projected within the central regions of 
NGC 4649, where the chance of finding a background source is low (due to the small solid 
angle involved). One might be concerned that the GC optical identification of one of the hard 
sources might actually be a background AGN. However, the identified GC has very normal 
red optical colors, which would be atypical (but not unheard of) for an AGN. This suggests 
that a small portion of the GC and non-GC X-ray sources may have strongly self-absorbed 
spectra. 

In some early- type galaxies, "supersoft" sources have been found which have no X-ray 
emission at energies above 1 keV (e.g., Sarazin et al. 2000). These sources can be fit with 
~75 eV blackbody spectra. In terms of their spectra, these sources are similar to the many 
supersoft sources observed in our Galaxy and in M31 (e.g., Kahabka & van den Heuvel 1997). 
In general, the Galactic supersoft sources arc believed to be accreting white dwarf stars in 
binaries, with the luminosity due to steady nuclear burning. However, the supersoft sources 
found in distant early-type galaxies have very large bolometric luminosities which exceed 
the Eddington luminosity for a Chandrasekhar mass white dwarf. A similar bright supersoft 
source has been found in the bulge of M81 (Swartz et al. 2002). One hypothesis is that these 
sources contain intermediate mass (~ 10^ — 10^ M©) accreting black holes (Swartz et al. 
2002). It is unfortunate that there are no supersoft sources in the regions of these galaxies 
which are covered by the GC surveys, particularly as it has been suggested that intermediate 
mass black holes are produced in GCs (Miller & Hamilton 2002). 



5.3. X-ray Variability 

During these Chandra observations lasting typically ~10 hours, about 10% of the X- 
ray sources are found to vary significantly. Due to the weakness of the sources and some 
features of the detector used, only secular variations over the observation period (source 
turns on, turns off, or increases or decreases by a significant fraction) are detectable. Of the 
111 LMXBs in our full sample, 9 varied significantly, of which 4 were in GCs. In the Lx 
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complete sample, 4 sources varied, of which 2 were in GCs. In the entire full sample, 44% 
of the LMXBs are in GCs. Within these very limited statistics, there is no evidence for a 
strong difference in the variability of GC and non-GC X-ray sources. 

5.4. X-ray/Optical Correlations 

Among the LMXBs which were located in GCs, we also searched for correlations between 
their X-ray properties (Lx, i?21°, if31°) and the optical properties of the GCs in which they 
are situated. For example. Figure 13 shows the logarithm of the X-ray luminosity (logLx) 
plotted against the optical color (V — I) of the GC for sources in the Lx complete sample. 
No significant correlations were found between any of the X-ray properties and any of the 
optical properties discussed here. The probabilities that the properties were uncorrclatcd 
were all 24% or greater for the various pairs of properties. One possible problem with these 
tests is the relatively small sample of only 49 LMXBs in GCs. 

6. Spatial Distributions 

As discussed in § 1, the available GC samples for these galaxies only cover a small 
fraction of the areas of the galaxies. As a result of this spatial incompleteness and possible 
position-dependent selection effects for both the GC and LMXBs, we will not analyze the 
spatial distribution functions of the GCs or LMXBs. Instead, we will only address the 
relative distributions of GC and non-GC LMXBs, and of X-ray and non- X-ray GCs. Also, 
we will only consider the radial distribution of sources. Because the galaxies have differing 
eccentricities, it would be difficult to compare their azimuthal distributions. In order to 
improve the statistics by combining the data from all of the galaxies, we will scale the 
projected radius r of each source by the effective radius r^s for its host galaxy. The adopted 
values of Ves are given in Table 1. The references are given in the papers containing the 
X-ray data. 

6.1. Distribution of LMXBs 

We first compare the spatial distributions of GC and non-GC X-ray sources. Figure 14 
shows histograms of the spatial distributions of these two types of LMXBs for the Lx com- 
plete sample. There is a slight tendency for the GC LMXBs to be located at larger radii 
than the non-GC LMXBs, but it is not very significant. The median radii of the GC and 
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non-GC LMXBs are 0.65 and 0.49 Ves, respectively. Using the Wilcoxon or Mann- Whitney 
rank-sum test or the KS test, the probabihty that the two samples were drawn from the 
same radial distribution can only be rejected at the 93% level. In the full sample, this trend 
is even weaker. The individual galaxies in our sample have too few sources to allow useful 
comparisons. Kundu et al. (2002) found that the spatial distributions of GC and non-GC 
X-ray sources were similar in NGC 4472. 

6.2. Distribution of GCs 

Next, we compared the radial distributions of GCs either containing or not containing 
LMXBs. Figure 15 shows histograms of the spatial distribution of X-ray and non-X-ray 
GCs in the complete sample. There appears to be some tendency for the X-ray GCs 
to be located at smaller radii. The median radius of the X-ray GCs is 0.67 Ves, while that 
for the non- X-ray GCs is 1.00 Tefr. The rank-sum test indicates that the radii of the X-ray 
GCs are smaller, with a probability that they were drawn from the same distribution of 
0.3%. Similarly, the KS test finds that the probabihty that the distributions are identical 
is about 0.5%. Similar results are found for the full sample, and in each of the individual 
galaxies, although the statistics for individual galaxies arc poor. Figure 16 shows histograms 
of the spatial distribution of X-ray and non-X-ray GCs in NGC 4365 and in NGC 4649. In 
NGC 4649 and in NGC 4697, the ranges of radii covered (in terms of Vcs) are both small, 
in NGC 4649 because the value of r^s in large in angular units (Table 1) and in NGC 4697 
because the GCs are all in an annulus at larger radii (Figure 1). 

Thus, there is at least marginal evidence that the X-ray GCs are more centrally located 
than the non-X-ray GCs. Kundu et al. (2002) found a similar effect in NGC 4472. Since the 
X-ray globular clusters are brighter and redder than the non-X-ray GCs, this weaker radial 
trend might be due to these other stronger correlations, and to a radial trend in the absolute 
magnitudes or colors of GCs within our sample. Indeed, we find that there is a correlation 
between GC optical luminosity or color and radius, with the brighter and redder GCs being 
more likely to be found at smaller radii. The optical luminosity gradient is certainly, at 
least in part, a selection effect, since it is harder to sec faint GC near the center of the host 
galaxy. Other selection effects might be introduced by differing absolute magnitude limits 
and differing radial coverage in each galaxy. In light of various selection effects in the raw 
data (e.g., the radial variation of incompleteness in the detection limits of GCs, which differs 
from galaxy to galaxy), and the radial variation in the intrinsic properties of GCs, such as 
the more spatially concentrated nature of red clusters with respect to the blue ones seen in 
most galaxies (e.g., Geisler, Lee, & Kim 1996; Kundu et al. 1999), it is not clear from this 
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data set whether the distance from the center of a galaxy independently affects the efficiency 
of LMXB formation in a GC. Crude statistical tests with this sample suggest that the radial 
variations due to selection effects and intrinsic properties might indeed be able to reproduce 
the observed radial gradient. 

7. Conclusions 

We have studied the connection between LMXBs and GCs in a sample of early-type 
galaxies. Ultimately, it is hoped that such studies may provide information on the origin 
and evolution of LMXBs, and on the evolution of GCs. In general, the fraction of LMXBs 
associated with globular clusters (/x-gc ~ 10 — 70%) is much higher (by factors of 10^ — 10^) 
than the fraction of optical light in GCs. This result, which is also true for our Galaxy, indi- 
cates that stellar dynamical interactions in globular clusters are very effective at producing 
compact binary systems (Clark 1975; Fabian et al. 1975; Hills 1976). 

There is evidence that the fraction of LMXBs associated with globular clusters increases 
along the Hubble sequence from spiral bulges to SOs to Es to cDs. The increase is a factor of 
~6. For some time, it has been known that the specific frequency of globular clusters (the 
number per unit galaxy bulge optical luminosity) also increases along the same sequence 
(e.g., Harris 1991) by a similar factor. These two trends might suggest that there are two 
populations of LMXBs in early-type galactic systems: those formed in GCs, and those formed 
from field binary stars. The variation in the fraction of LMXBs in GCs might then result 
from the variation in the specific frequency of GCs. If the number of field LMXBs were 
proportional to the field optical luminosity in all galaxies, and the number of GC LMXBs 
were proportional to the number of GCs, one might expect the variation in the fraction 
of LMXBs in GCs with Hubble type to be considerably smaller than the increase in the 
specific frequency of GCs, which is not really the case. However, variations in the field 
stellar populations from galaxy to galaxy might result in variations in the number of field 
LMXBs per optical luminosity. We note that observations suggest that the total number 
and X-ray luminosity of LMXBs in galaxies increase in proportion to the total number of 
GCs rather than the optical luminosity (White et al. 2002). This suggests that most or 
all LMXBs are formed in GCs (Grindlay 1984; Sarazin et al. 2001; White et al. 2002), at 
least in the earlier type galaxies. Then, the variation in the fraction seen in GCs today 
with galaxy type requires that a higher fraction of LMXBs have escaped from GCs in later 
type galaxies (spiral bulges) than earlier galaxies (gEs and cDs). However, the fraction of 
LMXBs which would have escaped from GCs in later type galaxies would be rather large, 
and might be inconsistent with models for the destruction of GCs or the stellar dynamical 
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escape of LMXBs. Future tests with larger samples may allow us to determine what fraction 
of LMXBs are made in GCs versus the field, and might provide independent information on 
the rate of destruction of GCs in galaxies. 

The fraction of globular clusters which contain X-ray sources appears to be roughly 
constant at /gc-x ~ 4% in different galaxy types (Kundu et al. 2002). Because of possible 
incompleteness of the GC samples at the faint end and because the probability of a GC 
containing a LMXB correlates with its optical luminosity, it may be more useful to give the 
number of high luminosity LMXBs per optical luminosity of GCs. In our complete sample, 
this value is about 1.5 x IQ-^ LMXBs per Lqj for Lx>lx 10^^ erg 8"^ (0.3-10 keV). The 
value rises to about 2.0 X 10"^ LMXBs per Lqj at lower X-ray luminosities [Lx ^ 3 x 10*^^ 
erg s"^). The value for the Milky Way is also very similar to this. The rate of occurrence 
of bright LMXBs in GCs may provide information on the duty cycle of quiescence vs. X- 
ray brightness for GC LMXBs, complementing information from the detection of quiescent 
LMXBs in Galactic GCs (e.g., Grindlay et al. 2001). 

The strongest trend which we have found is a correlation between the probability that 
a GC has an X-ray source and its optical luminosity. However, this correlation is consistent 
with a constant probability of finding a LMXB per unit optical luminosity; that is, it seems 
to result primarily from the larger number of stars in optically luminous GCs. Thus, this 
correlation doesn't provide any direct information on the physics of the formation and evo- 
lution of LMXBs in GCs. It would be very useful to correlate the probability that a GC 
contains a LMXB with the stellar density in the cluster in order to test the role of stellar 
dynamical interactions. The stellar densities could be estimated from the optical magnitudes 
and sizes determined from ffST observations. 

We find that LMXBs are more likely to be found in redder GCs, and that this trend is 
independent of the optical luminosity correlation. This seems to indicate that the evolution 
of X-ray binaries in GCs is affected either by the metallicity or age of the GC, with younger 
and/or more metal-rich GCs having more LMXBs. 

Models for the evolution of LMXBs suggest that the population should vary with time 
after the formation of the stellar population (White & Ghosh 1998; Wu 2001). White & 
Ghosh (1998) argued that the number of LMXBs would initially rise with time for ~ 1 
Gyr, which is set by the time scale for the secondaries in LMXBs to evolve. After this, the 
number of LMXBs should dechne slowly with time as the number of pre-LMXB systems 
decreases and the existing LMXBs evolve into millisec pulsar binaries. Given that the stellar 
populations in GCs are generally much older than 1 Gyr, we should be in the latter phase 
of evolution. Thus, a decline with age would not be unexpected. However, the arguments in 
White & Ghosh (1998) apply to LMXBs formed from the isolated evolution of a population 
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of initial binary stars; in GCs, new tight binary systems can be produced by stellar dynamical 
interactions. This might reduce (or even reverse) the rate of decline of LMXBs with age. 
More detailed models for the evolution of the population of LMXBs in GCs would be very 
useful. 

The number of LMXBs might also depend on the metallicity of the stellar population. 

A correlation between metallicity and the presence of LMXBs is also seen in GCs in our 
Galaxy and in M31 (e.g., Bellazzini et al. 1995). This might be due to the fact that higher 
metallicity stars are larger (Bellazzini ct al. 1995); however, the direct effect of this on the 
evolution of LMXBs docs not appear to be very strong (e.g., Ritter 1999). Larger stellar sizes 
might also promote tidal captures to form LMXBs (Bellazzini et al. 1995). Alternatively, 
metaUicity might affect the stellar initial mass function (Grindlay 1993) or have some other 
more complicated effect on the formation of LMXBs. 

There is no strong difference in the X-ray luminosities of GC and non-GC LMXBs in 
our sample. This disagrees with the result found by Angelini et al. (2001) for the galaxy 
NGC 1399, where the brightest LMXBs appeared to be positively correlated with GCs. On 
the other hand, our result is consistent with the analysis of the sources in NGC 4472 by 
Kundu et al. (2002). There is a weak tendency for the brightest LMXBs, whose luminosities 
exceed the Eddington luminosity for a 1.4 Mq neutron star, to avoid GCs. That may indicate 
that black hole X-ray binaries are somewhat less likely to be found in GCs, as seems to be 
true in our Galaxy. It may be difficult to retain massive BH binaries in GCs (Portegies Zwart 
& McMillan 2000). In any case, there arc some luminous LMXBs associated with GCs. If 
one assumes that these most luminous sources contain fairly massive BHs, it would appear 
that GCs can retain at least some massive BH binaries. 

There is no clear evidence that the X-ray spectra or variability of GC and non-GC X-ray 
sources differ. A much larger sample and repeated observations with Chandra would give 
stronger limits on the differences in the spectra or variability. 

We also find no evidence for a significant difference in the spatial distribution of GC 
and non-GC LMXBs. There is some evidence that X-ray GCs are more likely to be found 
near the centers of the galaxies than non-X-ray GCs, but this seems to be mainly due to 
an anticorrelation of optical luminosity with radius in our sample, plus the correlation of 
LMXBs with the more optically luminous GCs. Unfortunately, there are problems with the 
completeness of both the X-ray LMXB samples and optical GC samples. For the X-ray 
samples, the biggest problem is probably that most of the galaxies observed by Chandra so 
far have been X-ray bright galaxies, where diffuse emission from hot gas makes it difficult 
to detect fainter LMXBs near the centers of the galaxies. Hopefully, more observations of 
X-ray faint early-type galaxies will resolve this problem. Repeated Chandra observations 
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would also allow variability to be detected in more cases, and would allow fainter persistent 
sources to be detected in the combined exposure. The main limitation in the optical samples 
of GCs is that not all of the best galaxies have been studied with the HST, and the field 
of view of a single WFPC2 observation is smaller than the angular size of nearby early- 
type galaxies. Multiple observations with the WFPC2 or with the ACS may provide more 
complete samples of GCs for more galaxies. 
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This preprint was prepared with the AAS macros v5.0. 



Table 1. Sample Early- Type Galaxies 







Nh 


Tea 


D 


X-ray Exp. 


Count Rate Limit 


Lx Limit 


References 


Galaxy 


Type 


(10^0 cm-2) 


n 


(Mpc) 


(ks) 


(10-4 s-i) 


(10^^ erg s-i) 


LMXB 


GC 


NGC 1553 


SO 


1.50 


63 


18.5 


23.2 


4.3 


12.3 


BSI 


KWb 


NGC 4365 


E3 


1.63 


50 


20.4 


40.4 


2.7 


10.6 


SSI 


KWa 


NGC 4649 


E2 


2.13 


82 


16.8 


36.8 


3.1 


7.0 


RSI 


KWa 


NGC 4697 


E6 


2.14 


72 


11.7 


39.4 


2.6 


2.7 


SIB 


K 



References. — BSI (Blanton, Sarazin, 8z Irwin 2001); K (J. Kavelaars 2000, private communication); KWa (Kundu 
& Whitmore 2001a); KWb (Kundu & Whitmore 2001b); RSI (Randall et al. 2003); SIB (Sarazin et al. 2000, 2001); SSI 
(Sivakoff, Sarazin, & Irwin 2003) 
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Table 2. X-rays Sources in GC Survey Regions 





Src. 


Count Rate 


Lx (0.3-10 keV) 










Ml 


V -I 


Lx 


Galaxy 


No. 


(10-4 s-i) 


(10^'^ erg s-i) 


H2l0 


H31° 


Var.? 


GC? 


(mag) 


(mag) 


Comple 


NGC 1553 


3 


15.07±2.67 


43.3 


-0.34 


-0.35 












NGC 1553 


4 


7.64±1.94 


21.9 


-0.40 


-0.47 




V 


-10.63 


1.054 


V 


NGC 1553 


5 


9.04±2.08 


25.9 


-0.20 


-0.13 












NGC 1553 


7 


5.47±1.62 


15.7 


-0.66 


-0.11 












NGC 1553 


8 


8.46±2.03 


24.3 


-0.14 


-0.59 












NGC 1553 


11 


4.56±1.50 


13.1 


+0.09 


-0.69 












NGC 1553 


16 


40.04±4.22 


114.9 


-0.10 


-0.54 












NGC 1553 


17 


66.51±5.41 


190.9 


-0.15 


-0.65 


V 










NGC 1553 


18 


7.24±1.83 


20.8 


-0.09 


-0.51 












NGC 1553 


19 


4.46±1.52 


12.8 


-0.57 


-0.67 












NGC 1553 


23 


6.14±1.68 


17.6 


-0.43 


-0.72 






-12.17 


1.088 




NGC 4365 


1 


37.34±3.18 


152.5 


-0.17 


-0.53 












NGC 4365 


2 


14.23±1.98 


58.1 


-0.30 


-0.34 












NGC 4365 


3 


48.89±3.61 


199.6 


-0.13 


-0.40 












NGC 4365 


4 


8.86±1.56 


36.2 


+0.18 


+0.06 












NGC 4365 


5 


4.59±1.16 


18.7 


-0.42 


-0.42 






-9.26 


0.936 


V 


NGC 4365 


6 


5.94±1.35 


24.3 


+0.67 


-0.77 












NGC 4365 


7 


13.66±1.97 


55.8 


-0.26 


-0.23 












NGC 4365 


8 


3.37±1.01 


13.8 


-0.13 


-0.28 






-11.18 


1.069 


V 


NGC 4365 


9 


7.11±1.43 


29.0 


+0.31 


-0.11 




V 


-9.24 


1.038 


V 


NGC 4365 


10 


5.40±1.28 


22.0 


+0.44 


-0.06 






-8.51 


1.056 




NGC 4365 


11 


3.53±1.04 


14.4 


+0.01 


-0.41 












NGC 4365 


12 


7.18±1.45 


29.3 


+0.22 


-0.25 






-10.14 


1.241 




NGC 4365 


13 


4.16±1.14 


17.0 


-0.76 


-0.43 




%/ 


-9.95 


0.926 




NGC 4365 


16 


5.74±1.30 


23.5 


-0.11 


-0.42 


V 










NGC 4365 


17 


3.22±0.97 


13.1 


+0.58 


-0.26 






-10.68 


1.071 


V 


NGC 4365 


18 


11.78±1.80 


48.1 


+0.05 


-0.28 






-11.37 


1.125 




NGC 4365 


22 


3.93±1.06 


16.0 


+0.24 


+0.02 












NGC 4365 


25 


6.86±1.38 


28.0 


+0.20 


-0.06 






-11.18 


1.057 




NGC 4365 


26 


2.69±0.89 


11.0 


+0.33 


+0.00 






-8.64 


1.226 




NGC 4365 


29 


2.65±0.89 


10.8 


+0.24 


-0.01 


V 










NGC 4365 


30 


6.67±1.36 


27.2 


-0.37 


-0.49 












NGC 4365 


31 


2.66±0.89 


10.9 


-0.05 


-0.88 


V 










NGC 4365 


32 


4.05±1.08 


16.5 


+0.26 


+0.02 












NGC 4365 


33 


5.03±1.51 


20.5 


+0.41 


-0.76 












NGC 4365 


34 


2.94±0.93 


12.0 


-0.59 


-0.71 






-10.00 


1.205 




NGC 4365 


37 


2.79±0.92 


11.4 


-0.47 


-0.48 






-10.13 


1.290 




NGC 4365 


38 


18.65±2.27 


76.1 


-0.38 


-0.92 












NGC 4365 


40 


3.60±1.03 


14.7 


-0.07 


-0.85 






-11. .38 


1.101 




NGC 4365 


42 


11.86±1.79 


48.4 


+0.08 


-0.22 












NGC 4365 


43 


3.95±1.06 


16.1 


+0.46 


+0.18 










V 


NGC 4365 


44 


3.22±1.07 


13.1 


-0.05 


-0.70 






-11.15 


1.172 




NGC 4365 


45 


4.21±1.19 


17.2 


+0.19 


+0.09 






-10.98 


0.958 




NGC 4365 


46 


2.62±0.86 


10.7 


-0.05 


-0.70 












NGC 4365 


47 


10.72±1.68 


43.8 


+0.04 


-0.63 










V 


NGC 4365 


48 


2.70±0.89 


11.0 


+0.15 


+0.00 


V 




-9.86 


1.180 




NGC 4365 


50 


6.29±1.31 


25.7 


-0.26 


-0.72 






-9.69 


1.120 





-28- 



Table 2 — Continued 





Src. 


Count Rate 


Lx (0.3-10 keV) 








Ml 


V — I 


Lx 


Galaxy 


No. 




(10^''' erg s~^) 


H21° 


H31° Var.? 


GC? 






Comple 


NGC 4365 


53 


5.85±1.26 


23.9 


+0.55 


+0.44 


V 


-10.32 


1.027 


V 


NGC 4649 


2 


33.52±3.46 


75.4 


-0.61 


-0.90 










NGC 4649 


3 


20.43±2.90 


45.9 


-0.32 


-0.47 


V 


-9.93 


1.081 




NGC 4649 


4 


20.80±2.92 


46.8 


-0.44 


-0.40 










NGC 4649 


5 


13.11±2.50 


29.5 


-0.22 


-0.68 










NGC 4649 


6 


21.49±3.01 


48.3 


+0.13 


-0.55 


V 


-11.19 


0.999 


V 


NGC 4649 


9 


15.94±2.57 


35.8 


-0.21 


-0.42 










NGC 4649 


12 


9.94±1.98 


22.4 


+0.09 


-0.23 










NGC 4649 


16 


4.06±1.28 


9.1 


-1.00 


-0.28 










NGC 4649 


18 


20.77±2.80 


46.7 


-0.06 


-0.23 


V 


-11.03 


1.116 




NGC 4649 


19 


9.43±2.00 


21.2 


-0.11 


-1.00 








^/ 


NGC 4649 


20 


15.20±2.36 


34.2 


-0.08 


-0.22 


V 


-10.22 


1.225 


V 


NGC 4649 


21 


6.38±1.69 


14.3 


-0.15 


-0.34 








V 


NGC 4649 


22 


5.15±1.43 


11.6 


+0.01 


-0.56 


V 


-11.23 


1.246 




NGC 4649 


23 


9.29±1.85 


20.9 


+1.00 


+1.00 


V 


-11.27 


1.149 




NGC 4649 


25 


16.21±2.36 


36.5 


+0.20 


-0.10 










NGC 4649 


26 


6.69±1.66 


15.1 


+0.21 


+0.03 








^/ 


NGC 4649 


27 


7.68±1.75 


17.3 


-0.14 


-0.57 








%/ 


NGC 4649 


28 


7.39±1.72 


16.6 


-0.44 


-0.11 


V 


-10.19 


1.085 




NGC 4649 


29 


11.72±2.10 


26.4 


+0.05 


-0.68 


V 


-10.28 


1.228 




NGC 4649 


31 


18.10±2.49 


40.7 


+0.00 


-0.55 


V 


-10.57 


1.220 




NGC 4649 


32 


7.53±1.77 


16.9 


-0.50 


-0.54 


V 


-10.41 


1.031 




NGC 4649 


33 


6.16±1.60 


13.8 


-0.42 


-0.82 


V 


-8.46 


1.335 




NGC 4649 


40 


5.23±1.45 


11.8 


+0.09 


-0.01 










NGC 4649 


41 


5.92±1.48 


13.3 


+1.00 


+1.00 








V 


NGC 4649 


42 


11.38±2.01 


25.6 


+0.09 


-0.40 


V 


-10.19 


1.217 




NGC 4649 


45 


6.42±1.52 


14.4 


+0.04 


-0.43 


V 


-8.95 


1.218 


V 


NGC 4649 


46 


6.96±1.63 


15.7 


+0.54 


+0.19 


V 


-9.50 


1.185 




NGC 4649 


47 


5.88±1.51 


13.2 


-0.70 


-0.54 


V 


-10.19 


1.282 




NGC 4649 


48 


8.74±1.72 


19.7 


+0.03 


-0.04 


V 


-10.60 


1.218 




NGC 4649 


49 


4.89±1.36 


11.0 


+0.44 


-0.71 










NGC 4649 


53 


10.12±1.94 


22.8 


-0.04 


-0.06 


V 


-10.21 


1.320 




NGC 4649 


57 


3.89±1.16 


8.7 


+0.03 


-0.53 


V 


-11.33 


1.173 




NGC 4649 


58 


9.33±1.76 


21.0 


-0.20 


-0.32 










NGC 4649 


61 


3.85±1.16 


8.6 


+0.31 


-0.86 










NGC 4649 


62 


4.41±1.21 


9.9 


+1.00 


+1.00 


V 


-9.77 


1.058 




NGC 4649 


63 


9.51±1.74 


21.4 


-0.48 


-0.47 


V 


-11.46 


1.209 




NGC 4649 


65 


3.75±1.18 


8.4 


-0.90 


-0.64 










NGC 4649 


68 


4.87±1.28 


10.9 


-0.28 


-0.44 










NGC 4649 


69 


38.18±3.38 


85.8 


-0.20 


-0.69 










NGC 4649 


70 


19.65±2.44 


44.2 


+0.11 


-0.50 


V 


-8.24 


1.227 




NGC 4649 


71 


4.21±1.25 


9.5 


+0.53 


+0.20 










NGC 4649 


72 


6.34±1.52 


14.3 


-0.29 


-0.38 










NGC 4649 


73 


3.35±1.11 


7.5 


+0.74 


+0.17 










NGC 4649 


74 


16.57±2.26 


37.3 


-0.09 


-0.18 










NGC 4649 


75 


41.29±3.53 


92.8 


-0.12 


-0.44 










NGC 4649 


80 


10.88±1.88 


24.5 


-0.01 


-0.58 


V 


-10.29 


0.964 
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Table 2 — Continued 





Src. 


Count Rate 


Lx (0.3-10 kcV) 










Ml 


V - I 


Lx 


Galaxy 


No. 


(10-4 s-i) 


(10^^ erg s-i) 


H2l0 


H3lO 


Var.? 


GC? 


(mag) 


(mag) 


Complete? 


NGC 4649 


85 


17.74±2.37 


39.9 


+0.01 


-0.59 










V 


NGC 4697 


53 


25.75±2.62 


26.8 


-0.17 


-0.54 




V 


-7.40 




V 


NGC 4697 


54 


3.61±1.10 


3.7 


-0.47 


-0.48 




V 


-11.23 






NGC 4697 


55 


2.70±0.88 


2.8 


+0.43 


+0.39 


V 


V 


-8.16 






NGC 4697 


57 


10.16±1.64 


10.6 


-0.04 


-0.34 












NGC 4697 


58 


4.21±1.08 


4.4 


+0.29 


-0.13 


V 










NGC 4697 


59 


3.87±1.02 


4.0 


-0.07 


-0.75 












NGC 4697 


60 


2.62±0.86 


2.7 


-0.19 


-0.13 












NGC 4697 


61 


13.68±1.89 


14.2 


-0.50 


-0.62 


V 


V 


-7.91 




V 


NGC 4697 


62 


8.56±1.64 


8.9 


-0.07 


-0.28 












NGC 4697 


63 


7.46±1.51 


7.7 


-0.44 


-0.26 












NGC 4697 


64 


17.34±2.24 


18.0 


-0.12 


-0.93 


V 


V 


-9.37 




V 


NGC 4697 


65 


26.98±2.68 


28.0 


+0.18 


-0.37 










V 


NGC 4697 


66 


7.86±1.44 


8.2 


+0.00 


-0.86 












NGC 4697 


67 


3.67±1.10 


3.8 


+0.03 


+0.32 




V 


-6.59 






NGC 4697 


68 


3.27±1.02 


3.4 


+0.23 


-1.00 












NGC 4697 


69 


23.51±2.59 


24.4 


-0.00 


-0.29 




V 


-7.50 




V 



Table 3. Source Statistics in Sample Galaxies 



Lx Complete Sample Full Sample 













/x-GC 


/gc-x 


NBoth/Li 






/x-GC 


/gc-x 


NBoth/Li 


Galaxy 


Type 




A^LMXB 


A^Both 


(%) 


(%) 




A^LMXB 


A^Both 


(%) 


(%) 




NGC 1553 


SO 


70 


11 


2 






1.0±0.7 


11 


2 


18t?2 


2.91^1 


1.0±0.7 


NGC 4365 


E3 


325 


30 


14 






1.9±0.5 


37 


18 




5-5_i'.3 


2.0±0.5 


NGC 4649 


E2 


445 


36 


19 






1.7±0.4 


47 


22 


47l^ 


4.9+1-; 


2.0±0.4 


NGC 4697 


E6 


263 


5 


4 


801^4 


1 5+°-'' 

-■-•"^-0.7 


1.6±0.8 


16 


7 


44+9 


27+0.9 


2.8±1.1 



I 

CO 
o 

I 
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Fig. 1. — Regions covered by the GC survey for each of the sample galaxies are supposed on 
Digital Sky Survey (DSS) images of the galaxies. The images are 8' x 8', which is roughly the 
field of view of the Chandra images used to detect sources. In each case, the inner circle has 
a radius of one effective radius. For NGC 1553, NGC 4365, and NGC 4649, the polygonal 
region is approximate area of the HST GC survey. In NGC 4649, the quarter-circle region 
to the northwest was excluded as it appeared to mainly contain objects associated with the 
companion spiral NGC 4647. For NGC 4697, the GC survey region is the annulus between 
the two outer circles. 
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1 — ' — ' — ' — I — 1 




M, (mag) M, (mag) 

Fig. 2. — Histograms of the number of globular clusters versus their absolute magnitude, 
Mj. The left panel is for the Lx complete sample, while the right panel is for the full sample. 
In each case, the upper histogram is for all of the GCs in the galaxies. The lower shaded 
histogram shows the GCs which contain identified LMXBs. The histogram bins are 0.1 mag 
wide. 




Fig. 3. — Histograms of the number of globular clusters versus their absolute magnitude for 
the galaxies NGC 4365 (left) and NGC 4649 (right). The notation is the same as Figure 2. 
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Fig. 4. — Cumulative distribution functions for the probability that a GC contains an X-ray 
source ("X-ray") in the Lx complete sample and for the optical luminosity of GCs ("Opt. 
Lum."). 
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100 




100 




Fig. 5. — Histograms of the number of globular clusters versus their optical color, V — I. 
The left panel is for the Lx complete sample, while the right panel is for the full sample. 
In each case, the upper histogram is for all of the GCs in the galaxies. The lower shaded 
histogram shows the GCs which contain identified LMXBs. The histogram bins are 0.1 mag 
wide. 




Fig. 6. — Histograms of the number of globular clusters versus their optical color for the 
galaxies NGC 4365 (left) and NGC 4649 (right). The notation is the same as Figure 5. 
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Fig. 7. — Scatter plot of the absolute magnitudes M/ and optical colors F — / for the GCs. 
The open squares are non-X-ray GCs, while the filled squares are X-ray GCs in the Lx 
complete sample. 
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Fig. 8. — Histograms of the number of X-ray sources versus their X-ray luminosity (unab- 
sorbed 0.3-10 keV). The upper histogram is for all of the smirccs in our full sample. The 
lower shaded histogram shows the X-ray sources identified with GCs. The histogram bins 
are 0.1 dex wide. The arrow marks the limiting X-ray luminosity above which the sample is 
complete. 




Fig. 9. — Histograms of the number of X-ray sources versus their X-ray luminosity (unab- 
sorbed 0.3-10 keV) for the galaxies NGC 4365 (left) and NGC 4649 (right). The notation is 
the same as Figure 8. 
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Fig. 10. — Plot of the Galactic absorption-corrected X-ray hardness ratios, H21° and H31^, 
for the X-ray sources with at least 20 net counts in the Lx complete sample. The open 
squares are non-GC sources, while the filled squares are sources identified with GCs. The 
solid line and triangles show the hardness ratios for power-law spectral models, where the 
triangles indicate values of the power-law photon number index of F = (upper right) to 
3.2 (lower left) in increments of 0.4. The error bars at the upper left give the approximate 
uncertainties for an average source (one with about 50 net counts). 
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Fig. 11. — Plots of the X-ray hardness ratios if21° and if31° for the X-ray sources in the 
galaxies NGC 4365 (left) and NGC 4649 (right). The notation is the same as Figure 10. 
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Fig. 12. — Histograms of the number of X-ray sources versus their Galactic absorption- 
corrected X-ray hardness ratios H21^ (right panel) and H31^ (left panel). Only sources with 
at least 20 net counts in the Lx complete sample are included. In both panels, the upper 
histogram is for all of the sources. The lower shaded histogram shows the X-ray sources 
identified with GCs. The histogram bins are 0.2 wide. The horizontal bar at the top of each 
figures shows the typical uncertainty for that hardness ratio (the error bar for a source with 
about 50 net counts). 
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Fig. 13. — X-ray luminosity of LMXBs located in GCs versus the optical color of the GC. 
All sources are in the Lx complete sample. No significant correlation is evident. 
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Fig. 14. — Histograms of the number of X-ray sources versus their projected radius from the 
center of the galaxy in units of the effective radius of the galaxy, reff. The upper histogram 
is for all of the sources in the Lx complete sample. The lower shaded histogram shows the 
X-ray sources identified with GCs. The histogram bins are A(r/refr) = 0.2 wide. 
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Fig. 15. — Histograms of the number of GCs versus their projected radius from the center 
of the galaxy in units of the effective radius of the galaxy, Ves- The upper histogram is for 
all of the GCs in our sample. The lower shaded histogram shows the GCs identified with 
X-ray sources in the Lx complete sample. The histogram bins are A{r/res) — 0.05 wide. 



-44- 




Fig. 16. — Histograms of the number of GCs versus their projected radius from the center 
of the galaxy in units of the effective radius of the galaxy, r^s, for the galaxies NGC 4365 
(left) and NGC 4649 (right). The notation is the same as Figure 15. 



